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We find examples of duality among quantum theories that are related to arithmetic 
functions by identifying distinct Hamiltonians that have identical partition functions at 
suitably related coupling constants or temperatures. We are led to this after first developing 
the notion of partial supersymmetry — in which some, but not all, of the operators of a 
theory have superpartners — and using it to construct fermionic and parafermionic thermal 
partition functions, and to derive some number theoretic identities. In the process, we also 
find a bosonic analogue of the Witten index, and use this, too, to obtain some number 
theoretic results related to the Riemann zeta function. 
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1. Introduction 

Arithmetic quantum theories — theories for which the partition function is related 
to the Riemann zeta function or other Dirichlet series - - have yielded compelling 
results in mathematical physics J[J@. In addition to providing physical derivations and 
interpretations of various number theoretic results, such as an understanding of the Mobius 
inversion formula based on supersymmetry, these theories possess, as does string theory, a 
Hagedorn temperature, and so these theories may serve as toy models for certain aspects 
of string theory. 

In this paper, I find further evidence of the usefulness of arithmetic quantum thories 
in these areas. Most importantly, I find examples of duality in arithmetic theories, that is, 
of distinct arithmetic theories at different couplings that in fact have the same partition 
function. Such dualities are analogous to the coupling constant dualities of string theory.0 
In addition, I find physical underpinnings for additional number theoretic results. The 
most interesting examples here rest on "partial supersymmetry," a property identified in 
this paper, in which some, but not all, of the creation and annihilation operators have 
superpartners. This leads not only to a physical basis for some number theoretic results, 
but also to a formulation of parafermions of non-integer order and to the identification of 
an analogue for the Witten index based on a bosonic grading. 



2. Review of Notation and Ideas 

We here review the notation and way of thinking that allow a connection to be drawn 
between arithmetic number theory and quantum field theory. 

Consider some non-interacting quantum field theories defined in terms of their creation 
and annihilation operators. Let b\ and bk denote bosonic creation and annihilation 
operators, and let fl and fk denote fermionic creation and annihilation operators. In 
addition, we label the prime numbers pu in order of increasing magnitude, so that pi = 2, 
P2 = 3, p 3 = 5, etc. 

For use throughout this paper, we define the bosonic Hamiltonian 

oo 

H B =uJ2^g(Pk)blb k , (2.1) 

k=i 

1 See Section VII for a fuller discussion of this comparison. 
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and the fermionic Hamiltonian 

oo 

H F =u;J2 l °g(Pk)flfk ■ (2.2) 



fc=i 



We label the states of such theories using a Godel numbering. Let us denote the vacuum 
state |1). For the bosonic theory, we label the states! 

(6l 1 ) n (6D ra "-( 6 L) r "l 1 > = l^> > ( 2 - 3 ) 

where N = (pi) ri (P2Y 2 ' ' ' (PnY n ■ We may use a similar labeling for the states of the 
fermionic theory. Note that with this labeling, the Fock space of states for the bosonic 
field theory is labeled by the positive integers, while the Fock space of states for the 
fermionic field theory is labeled by the squarefree positive integers. One sees [0J that 

00 1 

PHB = E ^ = W (2-4) 



tre 

m 

m=l 

and 



*(-!)'.-"*= E 4£ ■ (2-5) 

m=l 

where n(m) is the Mobius inversion function 0. Supersymmetry tells us, via the Witten 
index [(§, that tr[(-l) F e~ f3( - HB+HF ^] = 1, and so 

V ^ = tr[(-l) F e"^] = - 1 — . (2.6) 



m=l 



3. Thermal Partition Functions 

In the preceding work, we have identified the bosonic thermal partition function with 
the Riemann zeta function, but only the graded fermionic partition function with a number 
theoretic function. On the other hand, the physics of a theory is governed by the genuine 
thermal partition function, and so to understand the thermodynamics of the fermionic 
field theory of Hp, we need to evaluate the partition function 

Z F {p y u) = tr e~ pHF . (3.1) 

2 For simplicity, we use the atypical convention that the creation and annihilation operators 
produce normalized states. 
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To calculate this partition function, we will introduce a notion of partial supersymmetry. 
This will then give us a physical and algebraic understanding for a particular zeta function 
identity.i 

Our approach is to start with the bosonic partition function and to cancel out just 
the right states so that we are left with the fermionic partition function. Suppose, then, 
that we define operators q% = (bk) 2 and ql = {b\) 2 , and in addition we define operators Ck 
and c\ which have the same effect as bk and b\ (respectively), except that Ck and c\ both 
square to zero. It is easy to see, then, that we can write the bosonic Hamiltonian 

H B = u;^log(p fc )4cfc + uj^22\og(pk)qlqk ■ (3.2) 

k k 

The two portions of this decomposition of Hb commute, and so 

tTe -PH B = tr [ e -/3« £ fc log( Pfe )ct Cfe] tr [ e -0« £ fc 2 log(p fe ) 9 t qk ^ ^ 

Because the Ck and c\ square to zero, the first term in this product is exactly the fermionic 
partition function. 

Notice that the q\, and q^ operators just appear as ordinary creation and annihilation 
operators, albeit associated with energy 21og(p/ £ ). (Due to our non-conventional 
normalization convention, there is no problem with how the resulting states are 
normalized.) Adding to the original Hamiltonian the expression ujJ2T=i^^ S(Pk)ftfk^ 
we are in a position to cancel out all but the and cl terms by taking a graded sum, by 
means of a Witten index type of result. It is this situation, in which we have fermionic 
superpartners for some of the bosonic creation and annihilation operators (the q\ and q^j 
but not others (the cl and Cfc) that we term partial supersymmetry. We thus find that 

Z f (J3,u) = tr[ e -^£ fc log(pfc)c fc Cfc ] 

= tr[(-l) F e"^^fe (log(pfe)c fe Cfe+21og(pfc)9 fe 9fe+21og(pfe)/ fc /fe:) ] (3-4) 
= tr[(-l) F e-^ HB+2(3H ^} . 

3 This partition function, like the parafermionic one of the next section, was obtained through 
number theoretic arguments in Thus it is our method in this section that is new, not the 
result. 
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The bosonic and fermionic pieces in the final expression in (|3.4j ) factorize into pieces that 
are easy to evaluate (see 0), leading to the number theoretic result 

Z F {M - ^ - ^ , (3.5) 

where s = /3u. Notice that we have obtained this number theoretic identity using a partial 
supersymmetry of a quantum mechanical theory, rather than the conventional explicit 
calculations of number theory using infinite products. 

Our use of partial supersymmetry to reduce a bosonic to a fermionic theory appears 
calculationally much like a BRST || cancellation. However, there is no gauge symmetry 
that is enforcing these cancellations, and so we stick to the term "partial supersymmetry." 
If one could identify a gauge invariance underlying this result, then one could perhaps recast 
these results as proper BRST cancellations. Such a result would undoubtedly provide 
further deep insights into the structure of these number theoretic quantum models. 



4. Parafermions of Arbitrary Order 

A natural generalization of the above results is to consider 

Z r {(3,u)=ti[{-l) F e-^ HB+rHF) ] . (4.1) 

As we know, Z\ = 1, as the Hp terms cancel all the Hg terms, yielding the Witten index 
of theory, while, as we just saw, Z2 = Zp, because the fermionic terms cancel all quadratic 
and higher powers of the original bosonic creation operators. 

For higher integers r, we can identify a partial supersymmetry as well. Let us define 
x\ and Xk such that their r th powers, but no lower powers, vanish as operators, and let 
r l = (^!) r an d r k =b r k . Thus 



H B = uJ^2(log(p k )xlxk + log(p r k )r{r k ) . (4.2) 



fc=i 



Then the fermionic terms in the expression for Z r (O) cancel the r\ and r k terms in Hb 



k 

and so 

Z r (J3, u) = tr[(-l) F e - p{HB+rHF) ] = tve- 0Hr , (4.3) 

where H r = u X^fcLi ^°E.(Pk)XkXk is the Hamiltonian for order r parafermions for the 
logarithmic spectrum we are studying. We can evaluate this expression, as the bosonic 
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and fermionic pieces of (|4.1| ) are non-interacting and hence factorize. One sees that the 
partition function for parafermions of order r, i.e., objects that are subject to an exclusion 
principle that states that no more than r — 1 parafermions can have the same quantum 
numbers, is given by 

What is compelling about this result is that we have found it using partial supersymmetry. 

One knows, of course, that the bosons are the infinite order limit of parafermions. 
It is interesting to see how that appears in this context. As we take r to infinity, the 
energies of states with fermionic excitations go to infinity, and so the fermionic creation and 
annihilation operators decouple 0, leaving behind the purely bosonic partition function. 
Thus the recovery of the bosonic theory at infinite r is, in our language, a decoupling effect. 

One of the most intriguing features of the partition function Z r is that it is well- 
defined for non-integer values of r. We do not simply mean that non-integer values of r 
can be plugged into the zeta function expression ( |4 . 4[ ) . Rather, we can put non-integer 



values of r into ( |4.1| ) quite sensibly and naturally — after all, why should the parameters of 
a Hamiltonian be integral — and this yields, from an expression with ordinary bosons and 
fermions, a partition function for exotic parafermions, parafermions of non-integral order. 
Thus, although there is no simple exclusion principle that one could state for such exotic 
parafermions, nor a simple Hamiltonian directly in terms of creation and annihilation 
operators for such exotic parafermions, we can obtain the thermal partition function quite 
easily from 

Z r ((3,u) =tv[(-l) F e-^ HB+rHF) ] . (4.5) 

What is the meaning of such theories? On the one hand, the partition function 
tells us what the theory is, so any quest for further interpretations might be viewed as 
superfluous. On the other hand, understanding what non-integral parafermions mean in 
a practical sense is physically intriguing. It is worth observing that the ability to vary 
the parameter r continuously in ( |4.3| ) means that we have a way to vary continuously 
from conventional parafermions of one integer order to another. This should enable one 
to connect results in the different conventional parafermion theories, presumably through 
analyticity or topological arguments. 

Note that for parafermions of order r real and less than 1, there is a temperature 
less than the Hagedorn temperature [|1] where the thermal partition function vanishes, 
indicating that the physical sense of these theories is limited to r > 1. This is no surprise. 



5 



At r = 1, the corresponding order 1 parafermions obey an exclusion principle such that 
states with one or more excitations are forbidden, which means that only the vacuum state 
is physical. (Indeed, r = 1 is simply the case of the Witten index calculation, in which all 
states but the vacuum state are canceled out.) Phrased in this way, parafermions of order 
r < 1 seem clearly to be physically untenable, as well as mathematically troublesome. 

5. Supersymmetry and Boson Number 

Recall from |ij that the mathematical interpretation of the supersymmetric Witten 
index result that 



IS 



= • (^) 

' m s Lis) 

while we also know from (|5 that the graded fermionic partition function is the inverse of 
the thermal bosonic partition function. Now that we have just calculated the thermal 
fermionic partition function, we can ask what its inverse is. We will find the result 
physically, and give its corresponding mathematical interpretation. 

Supersymmetry, we know, pairs bosonic and fermionic states. Consider the 
Hamiltonian Hs = Hb + Hp. Supersymmetry pairs each bk with and each bl with 
fl- Let us define (— 1) Nb to be the operator that has value +1 for Fock space states 
with an even number of bosonic creation operators, and —1 for Fock space states with an 
odd number of bosonic creation operators. Note that (— 1) Nb commutes with ft and ft- 
Because the supercharges of Hs involve terms like b\fk and f\bk, all states with non-zero 
energy come in superpartnerships with an equal number of even and odd Nb states. Only 
the zero energy states need not be paired in this way, as they are the supersymmetry 
singlet states. Consequently, for the theory with Hamiltonian Hs = Hb + Hp. we see 

tr[(-l) NB e- f3{HB+HF) } = 1 , (5.3) 

which is analogous to but distinct from the familiar Witten index result for tr[(— l) F e~ /3Hs ]. 
Since Hb and Hp commute with each other, we see that 

tre ~ l3HF = tr[(-l)"B e -/3HB] ' (5 ' 4) 
6 



The inverse of the thermal fermionic partition function is the graded bosonic partition 
function. 

The mathematical translation of this is as follows. Using the Godel numbering 
described above, we can represent (— 1) Nb by the arithmetic function ifj(m), where ifj(m) 
is totally multiplicative (i.e., , 0(m)'0(n) = ip{mn) for all natural numbers m and n) and 
with i^>(p r ) = (— l) r for any prime number p. Thus, using the connection with Dirichlet 
convolution describe in QTJ] , we see that 



where 1 is the identity function for Dirichlet convolution, l(m) = 5 m ,o, or, equivalently, 



Thus we have used our thinking about number theory to recognize that, generically, 
we can introduce a bosonic grading quite analogous to the fermionic grading, and this 
has led us physically to another number theoretic identity, namely ( |5.6|) . We note that 
the bosonically graded index is topological in nature, unaffected by deformations of the 
parameters of a theory as long as long as the supercharges remain odd under (— 1) Nb , and 
thus it is reliably calculable even by approximation methods. Supersymmetry is unbroken 
if this bosonic index is non-zero. Note, too, that (— 1) F corresponds to the arithmetic 
function /i, and (— 1) Nb to the arithmetic function One learns from supersymmetry 
that the inverses of these functions under Dirichlet convolution obey a nice symmetry: 
/U -1 = 1*01, while if) -1 = It would be interesting to find other pairs of functions with 
this type of partnering when taking the inverse, and to see if there is a supersymmetric 
reason underlying such a relationship. 

6. Graded Parafermion Partition Functions 

Our discussion of graded partition functions leads us naturally to consider graded 
parafermion partition functions. Here we present some results, but leave the detailed 
calculation to the readers, and reserve speculation on these results for elsewhere. 

Consider the parafermions of order r discussed previously. Define the operator fj such 
that a one-parafermion state has ^-eigenvalue 77, where rf = 1. Then 




(5.5) 




OO 



(5.6) 




(6.1) 
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using our representation of parafermions ( f4.3| ) and partial supersymmetry. Note that the 
bosonic creation operators also have ^-eigenvalue 77, while the fermionic creation operators 
are r)-neutral. 

There is an alternative graded partition function that one can naturally define. Let £ 
be an operator such that a one-parafermion state has ^-eigenvalue £ where £ r = —1. Then 
using our method of partial supersymmetry and the notion of the bosonic index ( |5.3| ), we 
find that 

ti[ie- fiHr ] = tr[|e-^ Hs ] tv[e- prHF ] . (6.2) 

The convention, again, is that the bosonic creation operators have ^-eigenvalue £, while 
the fermionic creation operators are ^-neutral. 

The 77- and ^-graded partition functions can be written in these factorized forms 
involving products of purely bosonic and purely fermionic terms because of existence of 
a partial supersymmetry in these expressions, which leads to cancellations along the lines 
of the Witten index in fl6.1fl and of the bosonically graded index (|5.3| ) in (|6.2j). Other 
gradings would not admit such a factorization. 

7. Duality 

The identities we have obtained above provide relationships for partition functions 
(thermal and graded) at different temperatures. This offers an interesting possibility, akin 
to duality in super symmetric and string theories: that we can identify two distinct theories 
at distinct temperatures that nonetheless give rise to the same thermodynamics. We will 
in fact realize this possibility, with an example in which the Euler totient function plays a 
central role. 

Before proceeding further, it is worth commenting on our use of the term "duality." In 
string theory, this term is most conventionally used to refer to relationships that directly 
rate the weak and strong coupling regimes (roughly speaking, replacing a coupling constant 
g by 1/g). Here, we use the term in a somewhat broader sense, to refer to identities 
among theories at different coupling constants in general. Although the relationships we 
identify do not simply involve inverting a coupling constant, they can be iterated to related 
weak and strong couplings. Furthermore, even without such iteration, the identities we 
find between arithmetic quantum theories can clearly play a role analogous to the role 
played by weak-strong coupling dualities in string theory, establishing the equivalence of 
superficially distinct theories, and thus allowing the use of one formulation of a theory to 
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analyze another. Hence we employ the term "duality" to characterize the identities among 
arithmetic quantum theories that we establish below. 

To begin our analysis, then, consider the thermal partition functions associated with 
the original Hamiltonians Hb and H F , which we will denote as Zb and Z F , respectively. 
Then ( |3.5|) can be interpreted as 

Z F (fl)Z B (2fl) = Z B (fl) . (7.1) 

Alternatively, if we denote the fermionic partition function graded with (— 1) F as A. F , we 
have 

A F (2fl) = Z F (fl)A F (fl) . (7.2) 

What is interesting about these expressions is that they relate partition functions at 
different temperatures. There are two interpretations possible here: 

(a) In (|7.1|) , we have a non-interacting mixture of two systems, each at its own 
temperature, equivalent to a third system at its own temperature; 

(b) In ( |7.2| ), we have a mixture of two systems thermalized with each other at one 
temperature equivalent to a different system at a different temperature. 

Note that in (b), not all the partition functions are thermal (some are fermionically 
graded), while in (a) we have an unthermalized mixture of two systems. Ideally, we would 
like a thermalized system at one temperature equivalent to a thermalized system at a 
different temperature. It turns out that this is possible to achieve. 

Consider the Euler totient function, defined by (p(m) = J2d-(d m)=i 1> i- e -' ^ ^ s the 
number of positive integers less than m that are relatively prime to m. Then it is a 
standard result in number theory that 



oo 



C(.-1) = CWU^ . (7.3) 

t J 't i l ' ' 



We can interpret this physically as follows. Let us imagine two theories: one is governed 
by the original Hamiltonian Hb from (|2.1|) above. The other has the same spectrum as 
Hb, but has degeneracies given by the totient function; it has <p(m) states with energy 
culog(m). We call the Hamiltonian for this theory H^ and its partition function 
Then we see that, in terms of the dimensionless parameter s = flu, 



Z B (s - 1) = Z B (s)Z ip (s) 
9 



(7.4) 



In other words, a thermalized mixture of two ideal gase, governed by the Hamiltonians Hb 
and H v , is equivalent simply to a higher temperature gas governed by Hb- Alternatively, 
we can phrase this as a coupling constant duality. The Hamiltonian Hb at coupling 
constant u — ^ and the Hamiltonian H^ at coupling constant g have the same partition 
function at temperature T = (3~ 1 . The presence of this duality indicates an important 
additional similarity connecting arithmetic quantum theories and string theories. 

The duality described above is clearly intriguing. There are two especially interesting 
directions to consider. 

(1) These expressions give a way to shift temperature via coupling new systems. Let 
us call a system governed by Hamiltonian Hb a Riemann gas, and one governed by H<f> 
a totient gas. Then what we have found is that if we start with a Riemann gas at one 
temperature, and then mix it with a totient gas at the same temperature, that this is 
equivalent simply to taking the Riemann gas and increasing its temperature a certain 
amount. Note that we can iterate this process by mixing in additional totient gases, and 
thereby shift the temperature further. This gives us an alternative to standard methods 
such as the renormalization group for thinking about how to shift temperatures; relating 
such distinct methods is potentially powerful. The most interesting outcome would be 
to determine a way to generalize this notion of obtaining a temperature shift by mixing 
in a new quantum system, finding ways to apply it to theories not so closely tied to the 
Riemann zeta function. 

(2) We have now identified, among the class of arithmetically-based quantum models, 
theories that exhibit duality behavior in their partition functions, with the best example 
being the duality (|7.4|) occurring for the theory with Hamiltonian Hb- We also know 
that this same theory exhibits a Hagedorn temperature. Both these traits — existence 
of duality and existence of a Hagedorn temperature — are distinguishing characteristics 
of string theory. This adds further support to the notion that the Hb theory and its 
relatives are related in a meaningful way to string theory, and thus that exploration of 
these arithmetic theories will prepare us to understand better the nature of string theory. 

We can expand on this idea a bit further. It is well-known that in string theory, the 
duality between large and small radius effectively means that there is a minimum distance 
in string theory. On the other hand, the Hagedorn temperature means there is a maximum 
temperature in a theory. Perhaps, in some similar fashion, the existence of a maximum 
temperature is associated with a temperature duality (recall that the duality we found 
could be understood as a temperature or as a coupling constant duality) in and of itself. 
These number theoretic field theories offer us an arena other than string theory in which 
to explore this question. 
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8. More with Partial Supersymmetry 

In earlier sections, we found that we could evaluate certain physically interesting 
partition functions, and simultaneously obtain some interesting number theoretic identities, 
using partial supersymmetry, in which some but not all of the bosonic creation operators 
had fermionic superpartners. There are other number theoretic results which appear as if 
they should have an interpretation via partial supersymmetry. We give here two particular 
examples to show the kinds of results that it appears ought to be possible to derive via 
partial supersymmetry, although to date no such derivation exists. 

As a first example, consider the function 

a a (m) = J2d a , (8.1) 

d\m 

where the summation is over all positive integers d that are factors of the positive integer 
m. The Dirichlet series of this function is related to the Riemann zeta function by 

C(s-a) = —-} S!±J. . 8.2 
v 1 C(s) ^ m s K ' 

3 v > 771 = 1 

Let us define H (Ta to be the Hamiltonian that has as its spectrum cr a (m) states of energy 
icTog(m). Then 

tr[e-^-S)^]=tr[(-l) F e-^-+^)] . (8.3) 

The right side involves a partial supersymmetric cancellation (note that none of the 
degeneracies go to zero) at temperature resulting in a theory a different temperature. 

A second representative case comes from re- writing the identity (|7.3|) of the preceding 
section as 

W) = =tr[e-0»-*>*»] tr[(-l)^e-^] . (8.4) 

In this context, the Euler totient function appears as the degeneracy factor resulting after a 
cancellation between bosons and fermions at different temperatures. This is a variation on 
the partial supersymmetry used above; here, if the bosons and fermions were at the same 
temperature, the cancellation would be complete. Mathematically, the right side of (|8.4| ) 
resembles a partition function for a supersymmetric Hamiltonian with non-supersymmetric 
boundary conditions. 

Needless to say, these two examples could both be written as identities involving a 
single temperature but different coupling constants. 
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The importance of highlighting these identities is to demonstrate what partial 
supersymmetry can offer. If these and a host of similar identities involving ratios of zeta 
functions and other Dirichlet series can be understood in this way, partial supersymmetry 
promises to be a valuable new tool for enriching our understanding of arithmetic number 
theory. 

9. Conclusions 

Arithmetic quantum theories have provided some of the most aesthetically pleasing 
examples of connections between mathematics and physics, and in this paper, we have 
found that these connections are even more compelling and worthy of study. 

Most importantly, we have found dualities among the partition functions for arithmetic 
quantum theories, dualities that can be interpreted as coupling constant or temperature 
dualities, demonstrating the equivalence of apparently distinct Hamiltonians. Our work 
has suggested that there may actually be reason to expect a connection between some types 
of dualities and the existence of a Hagedorn temperature, an area that warrants further 
investigation. With these dualities, it is possible to shift the temperature of a certain 
theory by mixing in additional quantum excitations at the original temperature; relating 
this interpretation to more familiar methods such as the renormalization groups should 
also provide useful insights into or constraints on the behavior of these theories. We are 
interested, too, in finding more useful applications of the bosonically graded topological 
index. 

On the technical side, the most significant result of this paper is the identification of 
the technique of partial supersymmetry. This has given a physical basis for a variety of 
number theoretic results, and we have identified further types of results that should yield 
to such an analysis. Indeed, one would expect this method to be useful beyond the realm 
of arithmetic quantum theories. 

This technique of partial supersymmetry has given us a physically-based formalism 
for interpolating from one Dirichlet series to another using ordinary bosons and fermions. 
Indeed, it is this idea that underlies our formulation of non-integer order parafermions. 
This is a potentially powerful idea, as it may allow the use of topological insights to 
study arithmetic number theory, by using continuity arguments to connect the behavior 
of different Dirichlet series. 
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It is our hope that the ideas introduced in this paper will be of value not only 
in understanding arithmetic quantum theories and directly related examples (physical 
ones such as string theory, mathematical ones such as Dirichlet series), but may also be 
extended to wider range of theories. There is no reason to think that the notion of partial 
supersymmetry or our particular approach to duality and its implications should not, in 
some fashion, be extended to relativistic, interacting theories with a variety of spectra 
(and not simply those theories so closely tied to the Riemann zeta function and other 
Dirichlet series), nor that we have exhausted the physical insights into number theories 
that consideration of quantum theories will provide. We hope further work will shed light 
on these possibilities. 

This research was supported in part by NSF Grant. No. PHY-9509991. 
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